Universal spectral statistics in quantum graphs 
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We prove that the spectrum of an individual chaotic quantum graph shows universal spectral 
correlations, as predicted by random-matrix theory. The stability of these correlations with regard 
to non-universal corrections is analyzed in terms of the linear operator governing the classical 
dynamics on the graph. 
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Fluctuations in spectra of individual complex quan- 
tum systems (e.g. classically chaotic systems) are uni- 
versal and can be described by the Gaussian ensem- 
bles of random-matrix theory (RMT). This statement, 
promoted to a conjecture by Bohigas, Giannoni and 
Schmit 0, has been empirically confirmed in numerous 
experimental and numerical analyses 0,0,0- However, 
never so far has it been possible to demonstrate analyt- 
ically that spectral fluctuations of individual 5] chaotic 
systems obey RMT statistics. Important progress has re- 
cently been made within the framework of periodic orbit 
theory: summing over orbit pairs of nearly identical ac- 
tion it became possible to prove universality (agreement 
with thepredictions of RMT) first for a few 0, Q and 
then all 9] coefficients of the short time expansion of the 
spectral form factor. Unfortunately it is not known how 
to advance this semiclassical approach into the regime of 
times larger than the Heisenberg time tn — (AE is 
the mean level spacing). 

Here, we apply different theoretical concepts to prove 
universality (including the long time regime) for a fam- 
ily of chaotic quantum systems, the so-called quantum 
graphs [loj . Quantum graphs differ from generic Hamil- 
tonian systems in that they are semiclassically exact (The 
density of states can be represented in terms of an exact 
semiclassical trace formula.), while they do not possess 
an underlying deterministic classical dynamics. Still they 
display much of the behavior of generic hyperbolic quan- 
tum systems; equally important, they are not quite as 
resistant to analytical approaches as these. 

In previous work, Berkolaiko et al. 0] developed a 
perturbative diagrammatic language to analyze the semi- 
classical periodic-orbit representation of spectral corre- 
lation functions beyond the leading ('diagonal') approx- 
imation. I n sp ite of the full knowledge of its building 
blocks [ill Il2l] a complete resummation of the pertur- 
bation series has so far been elusive (not to mention 
that such expansions are subject to the same limitations 
as the semiclassical approaches mentioned above.) In 
contrast, our present approach avoids diagrammatic re- 
summations. We rather build on two alternative pieces 



of input, both of which have been discussed separately 
before: i.) the exact equivalence of a spectral average 
for a quantum graph with incommensurate bond lengths 
to an average over a certain ensemble of unitary matri- 
ces 0,E3,[9, and ii.) the so-called color- flavor trans- 
formation which is an (equally exact) mapping of 
the phase-averaged spectral correlation function onto a 
variant of the supersymmetric cr-model. A subsequent 
stationary phase analysis then directly leads to the RMT 
correlation function corresponding to the symmetry of 
the graph. Finally, the spectrum of the 'massive' fluc- 
tuations around the saddle point contains quantitative 
information on the stability of RMT spectral statistics 
with regard to non-universal corrections. Deferring the 
discussion of other symmetry classes to a separate pub- 
lication |l(j| we will consider graphs which are invariant 
under both, time reversal and spin rotation. 

Let us begin by introducing our basic setting. A quan- 
tum graph consists of V vertices j connected by B bonds 
b. For the topology of the graph we assume that pairs of 
vertices are connected by at most one bond and that no 
bond starts and ends at the same vertex |l7|. We intro- 
duce 2B double indices (b, d), where d — 1,2 determines 
the (arbitrarily defined) direction of propagation along 
b = X,...,B. Boundary conditions on the graph are 
set by the fixed 2i3-dimensional unitary matrix Sbd.b'd' 
which describes the scattering of an incoming wave func- 
tion on bond b to an outgoing wave function on bond 
&'. Of course, Sbd,b'd' is non-vanishing only for bonds b 
and b' connecting at a common vertex j. Time-reversal 
invariance (T-invariance) implies that S T = af 11 Serf 11 , 
where af" — (<rf n )dd> are Pauli matrices in the space of 
directional indices. The complete dynamical information 
on the graph is carried by the 2B x 2B bond scattering 
matrix S(k) = T(k)ST(k). Here, the diagonal matri- 
ces T(k) contain the dynamical quantum phases picked 
up during propagation at fixed wave number k along the 
bonds: T{k)w,dd> = <W<W ex P where L b is the 
length of bond b and the two-fold replication in direc- 
tion space expresses the independence of the dynamical 
phases on the direction of propagation. The concise for- 



2 



mulation of this fact reads as T : S(k) — af 1T S T (k)af" . 
The prime signature of chaotic dynamics on the graph are 
strong non-Poissonian correlations in its discrete spec- 
trum {k n }. The latter is defined by the condition that 
S(k n ) has a unit eigenvalue or, equivalently, by the van- 
ishing of £(£;) = dct(l —S(k)) at k = k n . (This condition 
is equivalent 01 to the existence of an eigenvalue of the 
bond Schrodinger operator canonically associated to the 
scattering matrix. In this sense the spectrum {k n } is 
analogous to the discrete energy spectrum of a Hamilto- 
nian chaotic system.) Below, we will explore the two- 
point spectral correlation function 

R 2 (s) = A 2 (p(k + sA)p(k)) k -l 

where p{k) = ^2 n 5{k — k n ) is the spectral density and 

(. . . )k = liniA'^oc -ji $q dk(. . . ) an average over the wave 
number parameter k. (A = tt/BL the average level spac- 
ing and L = ~^2 h Lb/B the mean bond length.) 

Universal behavior of the spectral correlation func- 
tion (agreement with the prediction MT of RMT) 
can be expected if the corresponding classical system 
is chaotic (hyperbolic). What is the equivalent condi- 
tion on a quantum graph? An answer has been formu- 
lated by Tanner [3 (see also 0) in terms of the clas- 
sical probability Fu,b'd' = \Sbd,b>d> 2 = \S(k)bd,b'd' 2 to 
get from (b',d') to (b, d). This 'classical propagator' F 
has one eigenvalue Ai = 1, corresponding to equidistri- 
bution in bond space. The dynamics is mixing if, for 
large times, any initial probability distribution converges 
to this distribution, i.e. lim n ^ too (F n )bd,b'd' = 55 ■ This 
condition is met if all other eigenvalues |A2,...,b| < 1 lie 
inside the complex unit circle. However, mixing dynam- 
ics alone does not suffice to guarantee universality of a 
quantum grap h [l^. An additional condition proposed 
by Tanner 1 14| states that in the limit B — > oo, the spec- 
tral gap A g = max &S {2.... ! B}(l — |Af,|) be constant or, 
at least, vanish slowly enough A g ~ B~ a . Building on 
the so-called diagonal approximation (an approximation 
that obtains its asymptotics of i?2(p) for large values of 
p) Tanner conjectured that universal behavior should be 
expected for values of the gap exponent < a < 1. 

In the following we will show that quantum graphs 
indeed show RMT spectral correlations (provided the 
condition < a < 1/2 somewhat stronger than Tan- 
ner's is met.) We start out from the representation 
p(k) = A^ 1 - i jTmln£(fc+), k+ = k + iO of the den- 
sity of states in terms of the spectral determinant |Toj| . 
Using this formula, it is straightforward to verify that 
the two-point function assumes the form R 2 (p) = 

8^^f^Uo Re ^+'J-)>*' wh ere 



ak++p +b ) u(fc++P-b) 



(i) 



and P±h = (±s/2 - j±)A, P±f = (±s/2 + j±)A 



Our analysis will be based on the assumption that all 
bond lengths Lb are rationally independent. It has been 
shown [l3| that under this condition the average over the 
parameter k is strictly equivalent to an average over B 
independent phases e l_ 2~~ 1— > e ltph : 



(T[T(k)]) k = (^TO)]> , 



(2) 



where T is a smooth function of the bond-diagonal phase 
matrix T introduced above. 

To prepare the application of © to our present prob- 
lem, we represent the fraction of determinants in Q 
as a Gaussian integral, C — J d('ip,il')(^p(—S[ip,'^])) t p, 
where |20| 



S[ip,i/>]=ip- 



1 T(k) 
T(k) (ST+)t 



1p++1p- 



1 T(fe) + 
T(fc)+ ST- 



(3) 



Here, tp = {ip a ,s,x,d,b\ is a 16-B-dimensional supervec- 
tor, where a = ± distinguishes between the retarded 
and the advanced sector of the theory (determinants in- 
volving S and S' , resp.), s — f , b refers to complex 
commuting and anti-commuting components (determi- 
nants in the denominator and numerator, resp.), and 
x = 1,2 to the internal structure of the matrix kernel 
appearing in J2J). Defining a° as the Pauli-matrices in 
superspace, the matrices T± = T(2p±) i.e. T± are di- 
agonal matrices in superspace containing the bond ma- 
trices T(2p ±jb /f) in the boson-boson/fermion-fermion 
sector. Using that £(k + p) = det(l - ST(2p)) and 

det(l - ST(2p)) = det ( T ^ ^y^p)^ ) i one verifies that 
the Gaussian integration over all components of if) indeed 
yields the determinant 

As a second step, we subject the phase-averaged ip- 
functional to a duality transformation known as the 
color-flavor transformation |l5| . In a variant adapted to 
the present context [2l| , the transformation states that 



(exp(-S[&V])>* = (exp(-S'[* J *])) z , 



(4) 



where (•) = J dZ dZ sdct (1 — ZZ){-) and (matrix struc- 
ture in advanced/retarded space) 



S"[*,*] = ^ 



1 z 

Z T 1 



*1 + *2 



(ST+) + 
Z 



ST- 



*2- (5) 



Referring [2l| for a short discussion of the underlying 
technicalities, we here briefly explain the notation and 
the physical meaning of the transformation In 

©, *1,2 = {(4-1,2) a.s.t.b.d 

} are 165-dimensional inde- 
pendent supervectors, where the index, t — 1,2 ac- 
counts for the time-reversal symmetry of the model. 
Presently, all we need to know about the variables \E' 
and ^ is that they contain elements of tp and tp as 
their components, and fulfill 'J'i^ = ^2 T - Here, 
we introduced the fixed supermatrix r = af 11 ® To, 
where r = (E^a^ — iE^a^) (cf^ dl1 are Pauli matri- 
ces in the 'time-reversal' index t and direction index 
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d, E bh / S are projectors on the bosonic/fermionic sec- 
tor of the theory). The newly introduced integration 
variables, Z = bdiag [Z\, . . . , Zb) are 8-B-dimcnsional 
block-supermatrices with 8-dimensional entries Z^ = 
{Zb,ss',dd',w}- Finally, Z T = tZ t t~ 1 is, in a generalized 
way, transposed to Z , while Z and Z are independent. 





T 




T 






















Z 


























T 


- S 


T 





FIG. 1: On the physical interpretation of the color-flavor 
transformation. Explanation, see text. 

What is the physical significance of the transformation 
0? Fig. H shows a cartoon of the retarded (upper line) 
and advanced (lower line) wave function dynamics in the 
system. During propagation, both states pick up random 
scattering phases T (indicated by vertical dashed lines) 
and suffer scattering from one bond to the other (S- 
matrix). The rapid succession of these events implies wild 
fluctuations of the wave function amplitudes. Within 
the field theoretical context, this translates to uncontrol- 
lable fluctuations of the bilinears i'+ : s,x,d,be % '^ b '^+,s,x,d,b 
appearing in the action Q. In contrast, the field Z en- 
ters the theory as ~ * +, s ,t,d,bZb,ss> ,tt> ,dd'^ -,s',t' ,d',b, he. 
through structures that couple retarded and advanced 
field amplitudes (the 'vertical' ovals in the figure) . These 
amplitudes generally interfere to form slowly fluctuating 
entities (the basic principle behind the formation of uni- 
versal correlations.) This indicates that the Z-integral 
will be comparatively benign and, foreseeably, amenable 
to stationary phase approximation schemes. To pro- 
mote this expectation to a quantitative level, we inte- 
grate out the 'F's, thus arriving at the exact representa- 
tion (C) fc = J dZdZexp(-S[Z, Z]), 

S[Z, Z]=- str In (l - Zz) + -str In (1 - Z T Z) 

I (6) 
+ -str In (l - (5T_) t Z(ST+)Z r ) . 

As a first step towards a better understanding of the 
physics of this expression, let us consider its quadratic 
expansion, 

S (2) [Z,Z}=str {zZ -\z T Z - ^T^ZST+Z^ . (7) 

Void of non-linearities (terms of 0(Z 4 )), the action S^ 2 ' 
describes the un-interrupted propagation of two ampli- 
tudes along the same path in configuration space, i.e. the 
level of approximation underlying the diagonal approxi- 
mation in semiclassical periodic-orbit theory. 



This connection is made quantitative by noting that 
the action possesses a family of approximately (up to 
corrections of 0(B~ 1 )) 'massless' configurations, or 'zero 
modes' identified by SzS^ — SgS^ = 0. Upon substi- 
tution of the ansatz Zdd' — ddd'Zd (configurations not 
diagonal in direction space do not qualify as solutions, 
see the discussion of deviations below) these equations 
assume the form 

Z = Z T , (1-F)Z = 0, (8) 

where we have set the parameter matrices T± = l[22|. 
Owing to the fact that on a chaotic quantum graph 
the 'classical propagator' F has only one eigenvalue 
1, Eq. |JS} possesses the unique solution Zb,dd' — 
(2B)~ 1 f 2 5d,d'Y , proportional to the invariant equidistri- 
bution. We define Y T — tqY 7 "^ 1 , where the matrix 
To differs from r by the absence of the (now redundant) 
matrix af". Technically, the relation Y — Y T identi- 
fies (Y, Y) as generators of the orthosymplectic algebra 
osp(4|4). 

To explore the contribution of these modes to the the- 
ory, we set T±_b = 1 + iLbp± (provided the bond length 
fluctuations are not too large, Lb/L = O(B ), higher 
orders in the expansion may safely be neglected.), do 
the Gaussian integral over Y and differentiate w.r.t. the 
source parameters j± . As a result we obtain ^GOE,diag _ 
l/(7r 2 s 2 ) which agrees with the s ^> 1 asymptotics of the 
RMT correlation function. This is consistent with the 
observation that non-Gaussian contributions to the ex- 
pansion of S[Y] can no longer be neglected once s < 1. 

However, before going beyond the level of the Gaussian 
approximation, let us briefly consider the role of non-zero 
mode fluctuations. A glance at Eq. J2J shows that devi- 
ations from the first of the two equations in (JHJ are pe- 
nalized by a large action = 0(1). Upon integration, 
these modes produce a factor unity to the spectral deter- 
minant. Similarly, due to the absence of multiple connec- 
tivities and vertex loops, modes that are off-diagonal in 
the direction index d can be integrated out to give a factor 
unity. To explore the more interesting role played by de- 
viations from the equation (1 — F)Z = 0, let us expand a 
general configuration Z b ,dd> = <W Y^n=i Y mXm.bd in the 
basis of eigenf unctions Xm °f the operator F. Here, Y m 
are four-dimensional supermatrices obeying the symme- 
try Y m = Y^ and the identification Y\ = Y is understood. 
For T± = 1 the action of the 2B — 1 ~ B modes Y m> \ 

is given by |E m ( 1_ A m )strF m Y m > YE m str ^» y ™- 
Gaussian integration over these modes obtains a contri- 
bution ~ £ m (A/(l - A m )) 2 ~ 1/A 2 S to the spectral 
function. We conclude that the cumulative contribution 
of the 'massive' modes can safely be neglected provided 
£?A 2 -^ S ^s°° 23] . (The same holds true for higher or- 
der correlation functions, provided the order of the func- 
tion is smaller than J5.) 

Going beyond the level of the quadratic approximation, 
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we note that the saddle point equations SzS = S^S = 
of the full action ijBJ are still solved by the zero mode 
configurations While deviations from the zero modes 
continue to be negligible (as long as BA 2 g — > oo), the 
action of the latter now reads 

S[Y] = iBLstr (p_YY/(l - YY) - p+YY/(l - YY)j , 

where we have rescaled Y — > {2B) X / 2 Y . To represent this 
result in a perhaps more widely recognizable form, let 
us define the 8 x 8-matrix Q = ( \ \) ( \ _° x ) ( £ \)~ X . 
It is then straightforward to verify that the action S[Y] 
assumes the form of Efetov's^ action for the GOE cor- 
relation function 

1 7T 

S[Q] = ystr(Qe), (9) 

where e = — dia,g(p^ , p^ , p^ , p^ ) / A and the integration 
measure dQ — dYdY. 

Summarizing, we have proven Tanner's conjecture on 
universal spectral statistics in large chaotic quantum 
graphs. Our analysis was based on the assumption of (i) 
incommensurate bond lengths, (ii) the absence of multi- 
ple connectivities (iii) moderate bond length fluctuations, 
Lb/L < O(B), and (iv) weak scaling of the spectral gap 
A g ~ B a ,0 < a < 1/2 of the 'classical propagator' on 
the graph (which is stronger than Tanner's expectation 
< a < 1.) The conditions (i)-(iv) are met by various 
families of graphs 0, Q, 0] including the class of fully 
connected graphs with Neumann boundary conditions. 

We have enjoyed fruitful discussions with Fritz Haake, 
Sebastian Miiller, Stefan Heusler, and Peter Braun. This 
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Forschungsgemeinschaft . 



* gnutz@physik.fu-berlin.de 

* alexal@thp.uni-koeln.de 

[1] O. Bohigas, M.J. Giannoni, and C. Schmit, 

Phys. Rev. Lett. 52, 1 (1984). 
[2] T. Guhr, A. Muller-Groeling, and H.A. Weidenmuller, 

Phys. Rep. 299, 189 (1998). 
[3] F. Haake, Quantum Signatures of Chaos (2nd edition, 

Springer, Berlin, 2000). 
[4] For a discussion of some prominent exceptions to this 

rule, see J.P. Keating, Nonlinearity 4, 309 (1991), 

E.B. Bogomolny, B. Georgeot, M.J. Giannoni, and 

C. Schmit, Phys. Rep. 291, 220 (1997). 
[5] It has been known for some time |(| that ensembles of 

disordered systems exhibit RMT statistics upon disorder 

averaging. 

[6] K. Efetov, Super-symmetry in Disorder and Chaos, (Cam- 
bridge, 1997) 

[7] M. Sieber, and K. Richter, Phys. Scr. T 90, 128 (2001); 
M. Sieber, J. Phys. A 35, L613-L619 (2002). See, also 



I.L. Aleiner and A.I. Larkin, Phys. Rev. B 54, 14423 

(1996) . 

[8] S. Miiller, Eur. Phys. J. B 34, 305 (2003), D. Spehner, 
J. Phys. A 36, 7269 (2003), M. Turek and K. Richter, 
J. Phys. A 36, L455 (2003). 

[9] S. Miiller, S. Heusler, P. Braun, F. Haake, and A. Alt- 
land, to be published in Phys. Rev. Lett (ln9189). 
[10] T. Kottos and U. Smilansky, Phys. Rev. Lett 79, 4794 

(1997) ; Ann. Phys. 274, 76 (1999). 

[11] G. Berkolaiko, H. Schanz, and R.S. Whitney, 
Phys. Rev. Lett. 88, 104101 (2002); G. Berkolaiko, 
H. Schanz, R.S. Whitney, J. Phys. A 36, 8373 (2003); 
G. Berkolaiko, Waves in Random Media 14, S7 (2003). 

[12] S. Gnutzmann, B. Seif, Phys. Rev. E 69, 056220 (2004). 

[13] F. Barra and P. Gaspard, J. Stat. Phys. 101,283 (2000). 
The heuristics behind the proof is that for incommensu- 
rate bond length, k i— > { expi^r*, . . . , exp defines 
an ergodic flow on a B-torus. This means that the k- 
average can be replaced by a configurational average over 
the B-torus, i.e. an average over B independent phases. 

[14] G. Tanner, J. Phys. A 34, 8485 (2001). 

[15] M. Zirnbauer, J. Phys. A 29, 7113 (1996); M. Zirnbauer, 
in I.V. Lerner, J.P. Keating, and D.E. Khmelnitskii 
(edts.) Supersymmetry and Trace Formulae: Chaos and 
Disorder (Plenum, 1999). 

[16] S. Gnutzmann and A. Altland, to be published. 

[17] These technical assumptions simplify the analysis. We 
believe them to be physically irrelevant. 

[18] G. Berkolaiko, J. Phys. A 34, 319 (2001). 

[19] For an example of a mixing graph with non-universal 
spectral statistics see: G. Berkolaiko and J.P. Keating, 
J. Phys. A 32, 7827 (1999). 

[20] For a definition of the supersymmetric generalization of 
the determinant, 'sdet' and the trace, 'str', see |g. 

[21] Some details on the derivation of the action ©: start- 
ing from we introduce enlarged fields according to 
ip± -> *± = -^(ip±,ip±o-i ir ) T and i>± -> *± = 
-^{tp±,tp±(o'i 1T ® o"3 ))• The new fields depend on each 
other through *± = 9±t, where the orthosymplectic 
transposition r is defined in the main text. As a re- 
sult, the action assumes the form S = 2'5+iT^'+2 + 
2vj>_ 2 T t *_i + So, where So is the T-independent con- 
tribution and the indices 1,2 refer to the ^-components 
of 

To each of the B phase averages over e" l ' b € U(A^ C = 1) 
we now apply the color-flavor transformation in a variant 
(N c — l,Nf — 4), i.e. for one color and four 'flavors' 
(the latter labeled by t = 1, 2 and d — 1,2). Which in 
practice means flBl that the phase action gets replaced 
by S -> S' = 2*+iZ*_i+2f - 2 Z^ +2 + S where Zand 
Z are 8-B-dimensional block matrices containing 8x8- 
supermatrices Z b and Z h as independent sub-blocks. We 
finally use relations such as f _i2$+i = '5' +\Z T< S! _i to 
symmetrize the action and arrive at JKJ. 

[22] This simplification is justified because in practice (see 
below) T± = 1 + C(B- 1 ). 

[23] In deriving this result, we have neglected the coupling 
of the modes introduced by T±. However, for Lt/L = 
O(B ), this coupling leads to corrections to the action of 
0(B^ 1 ) which are negligible. 



